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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 2550 

DETERMINATION OF SHAPES OF BOATTAIL BODIES OF 
REVOLUTION FOR MINIMUM WAVE DRAG 
By Mac C. Adams 


SUMMARY 


By use of an approximate equation for the wave drag of slender 
Bodies of revolution in a supersonic flow field, the optimum shapes of 
certain hoattail bodies are determined for minimum wave drag. The prop- 
erties of three specific families of bodies are determined, the first 
family consisting of bodies having a given length and base area and a 
contour passing through a prescribed point between the nose and base, 
the second family having fixed length, base area, and maximum area, and 
the third family having given length, volume, and base area. The method 
presented is easily generalized to determine minimum -wave -drag profile 
shapes which have contours that must pass through any prescribed number 
of points. 

According to linearized theory, the optimum profiles are found to 
have infinite slope at the nose but zero radius of curvature so that the 
bodies appear to have pointed noses, a zero slope at the body base, and 
no variation of wave drag with Mach number. For those bodies having a 
specified intermediate diameter (that is, location and magnitude given), 
the maximum body diameter is shown to be larger, in general, than the 
specified diameter. It is also shown that, for bo.dies having a specified 
maximum diameter, the location of the maximum diameter is not arbitrary 
but is determined from the ratio of base diameter to maximum diameter. 


INTRODUCTION 


The wave drag of slender bodies of revolution having cross-sectional 
areas with a zero slope at the base was shown by Von Karman (reference l) 
to be given approximately by a double integral dependent only on the body 
shape and independent of Mach number. By use of Von Karman f s integral 
and the calculus of variations, several authors (references 1 to k) have 
treated the problem of determining optimum body shapes to give minimum 
wave drag. All these investigations have been concerned with either 
closed bodies or a body, such as a shell, having its maximum thickness 
at the base; however, none have treated bodies having boattails, and 
this problem is considered herein. 
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Ward (reference 5) has shown that bodies having a finite slope at 
the base give rise to drag terms in addition to Von Karman ! s integral 
and these additional terms include Mach number effects. The present 
paper shows , how ever, that the optimum bodies must have a zero slope at 
the base and consequently the additional terms vanish. The determination 
of the minimum-wave-drag bodies of revolution with boattailing then 
resolves itself to minimizing the same integral as used by the reference 
papers but with a more general treatment of the body profile. 

Although the analysis to follow is concerned with wave drag only, 
it should be remembered that the additional drag resulting from base 
pressure and skin friction is also of importance. A brief discussion 
of this additional drag is given in the concluding remarks. 


SYMBOLS 


D 


wave drag 


q 


dynamic pressure of free stream 


M 


free-stream Mach number 



Z 


body length 


distances made nondimens ional with respect to Z/2 and 
measured along body axis from midpoint of body 


r(x) 


body radius 


S(x) 


nondimens ional body cross-sectional area 



A 


specified cross-sectional area of body divided by (z/ 2)^ 


a 


body diameter corresponding to area A 


S, 


max 


maximum body cross-sectional area divided by (Z/2)^ 


B 


body- base area divided by (z/2) 2 
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D 


maximum body diameter 


drag coefficient 


W 1 / 2 ) S max/ 


V body volume divided by (Z/ 2)3 

c distance, divided by Z/2, from midpoint of body to loca- 

tion of specified area A 

e distance, divided by Z/2, from midpoint of body to loca- 

tion of S max 

\ Lagrange multiplier 


ANALYSIS 


(reference 5) has shown that 
a supersonic flow field is 


6|dx d| + 

2 / 

2iii f^S"U) 108,(1 - x)dx - i[s'(l)] 2 106e (1) 


where the body cross-sectional area and the coordinate distances have 
been written in nondimens ional form. The requirements in the derivation 
of equation (l) are that the body profile have no corners (a continuous 
slope) and that the rate of change of the cross-sectional area at the 
nose be zero (S'(-l) = 0). 

Optimum body having given length and base area and a contour passing 
through prescribed point between nose and base .- In certain practical 
problems it may be desired to have a minimum-wave-drag body of revolution 
which has a given length and base area and which has a contour that must 
pass through a specified diameter at a particular point between the nose 


By means of linearized theory, Ward 
the wave drag of a body of revolution in 
given approximately by 


/Z\2 


L 


-i 


S"(x)S"(|)log e |x - 


k 
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and base of the body. Such a situation would arise , for example, if a 
body of revolution were desired which would enclose some given rocket 
or jet engine. 

Although the analysis to follow determines the optimum body shape 
which is required to pass through three given points, that is, the nose, 
the base, and some in-between point, the manner in which the analysis 
could be generalized to determine optimum body shapes which must pass 
through any prescribed number of points will become apparent. 

Let it be supposed that the optimum body shape is given by the 
expressions 


S(x) = f(x) 
S(x) = g(x) 


(-1 = x = c) 
(c ^ X < 1) 


^ (2) 


Since the body must pass through three given points, the following condi- 
tions apply: 




f(-l) = 0 
f(c) = g(c) 
gU) = B 


= A > 
vj 


(3) 


The limitation of equation (l) allows no corners on the body profile; 
hence , 


f»(c) = g’(c) 


(4) 


and a further condition requires that f ? (-l) = 0. 
Next consider another body shape such that 


S(x) = f (x) + 6h(x) 
S(x) = g (x) + 7 k(x) 


(-1 = x ^ c) 
(c x % l) 


( 5 ) 
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where 


h(-l) = h(c) = k(c) = k(l) = h'(-l) = 0 
6 h*(c) = 7 k'(c) 


( 6 ) 


This second tody profile passes through the same three prescribed points 
as the optimum body profile and also satisfies the other conditions placed 
on the optimum body. However, since f(x) and g(x) define the optimum 
profile, the variations in the profile, 5 h(x) and 7 k(x), must neces- 
sarily produce an increased drag. It follows that h(x) and k(x) are 
completely arbitrary in form and that 5 and 7 are independent param- 
eters. Substitution of equations (5) into equation (l) gives the drag 
as a function of the two independent parameters S and 7. If the rear 
portion of the body is held fixed in its minimum-drag configuration, that 
is, 7=0, then 


5d(s,o) 

d5 


J 5=0 


= 0 



d| 



f"(x)loge|x 


l|dx + 



1 1 dx 


+ g’(D 



h" (x) logg (l 


x) dx 

(7) 


where 


h(-l) = h(c) = h * (- 1 ) = h'(c) = 0 

In a similar manner, the forward portion of the body can be held fixed 
in its optimum configuration, that is, 5=0, and 
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dD(0,7) 
by 


= 0 


7=0 


■ -L 1 k,,(5) 


d| 


f"(x)log e |x - | I dx + 


g"(x)log e |x - l| dx 


+ k«(l) 


f " (x ) log ( 1 - x) dx + 


-1 



g"(x)log e (l 


x) dx 


+ g ! (i) 



k"(x)log e (l 


x) dx - 


g'(l)k'(l)loge §^3 


( 8 ) 


where 


k(c) = k(l) = k'(c) = 0 


Partial integration of equations (7) and (8) yields 


and 



h(!) dS 




f»(x) 
x - I 


dx + 



= 0 


(9) 



k(ft) d! 

dl 2 


£ 8 (x) 


dx + 


-1 



+ 


g*(l) lim [k*(x)log e (l - x)| 

x->l 


g»(l)k‘(l)loge = 0 

2 V it 


(10) 
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Inasmuch as h(t) is an arbitrary function, equation (9) yields directly 


I dx + J dx = A l£ + A 2 (-1 = 1= c) (11) 


Equation (10) cannot be reduced so readily but requires more considera- 
tion. Since k(|) is also arbitrary, it can be chosen such that 
k ! (l) = 0 of sufficient order and then 


f dx + f dx = a 3 i + a 4 (c < 1 < 1) (12) 


Now suppose that k ! (l) is not zero but instead is some 
It then follows from equations (10) and (12) that g T (l) 
that is , the optimum body profile must have a zero slope 
It is interesting to observe that the Mach number effect 
disappears for the minimum-drag bodies. 

The simultaneous solution of integral equations (ll) and (12) will 
now give the optimum body profiles sought. These two equations are 
identical to those found in reference 2 in which minimum-drag closed 
bodies were considered. 

It is recognized that equations (ll) and (12) are analogous to the 
integral equations found in lifting-line theory. In the analogy f ! (x) 
and g’(x) correspond to the vortex distribution, and the right-hand 
sides of the equations represent the induced-normal- velocity distributions 
on the two-dimensional wings with no thickness. The cross-sectional area 
of the optimum body is therefore a linear combination of the surface 
potential distributions which give rise to the following normal- velocity 
distributions on two-dimensional wings: 


finite value, 
muet be zero; 
at its base, 
on the drag 


3 


NACA TN 2550 



where 


= A^c + A]^ C3 = A^c + A2 

C2 = A^ C|j. = -A-^ 


In addition, a circulatory solution is required in order that the body 
have a finite base area; that is, the analogous wing has a specified 
lift. 


Sketches (a) and (c) correspond to two-dimensional flat-plate wings 
with deflected flaps, and sketches (b) and (d) represent flaps rotating 
with a constant angular velocity. The surface -potential distributions 
for the four wings, as well as the circulatory solution, are all found 
in reference 6 and may be written 


*(a) 

*(c> 




(x 



cos"^c - 



cos 


-1 


(-C) 


■h (x - c)log e N 


(13a) 


(x - c)2loggN 

(I3h) 


(13c) 
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\jl - /l - x^ - (x - 2c) \/l - x^ coa"'*‘(-c) - (x - c)^loggN 

(13d) 

p 

9(r) = 2^ cos " 1 (- x ) (l3e ) 



where 


N = 


1 - cx 



and 9(p) represents the circulatory solution. The constants in equa- 
tions (13) are determined by applying the following conditions: 


S‘(-l) = S*(l) = 0 " 


S(c) = A \ 

S(l) = B 


m 


A further condition is needed to insure that the body have finite and 
continuous slope at the point C) that is, Cp = C3. 

The least-drag body, which has given length and base area and which 
has a contour that must pass through the point c where S = A, is then 
given by 


S(x) = 


— COP - "*"' 


A - — cos^C-c) 


(-c) 

\/l - x^ 

JTW 2 

(-C) 

Be 


(1 - cx) + ^ 


\l ~ 


1 - C c 


(x - c) + 


L (1 - =2)2 ' , (l . C 2)^J 


B 


(x - c) 2l oge N 4- - cos-l(-x) 


( 15 ) 
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Some typical profiles defined by equation ( 15 ) are shown in figure 1, 
where an arbitrary thickness ratio of approximately 20 percent has been 
chosen. All the profiles defined by equation (15) have infinite slope 
at the nose. However, the radii of curvature approach zero at the noses 
so that the noses appear to have finite angles. 

For vanishing base area, equation ( 15 ) reduces to the body profiles 
considered in reference 2. The optimum shell for given caliber and 
length is obtained for A — >B and c — >1. The result is 


S(x) 


= — x \fl - x 2 + — cos"-*-(-x) 

Jt V Jt V ' 


( 16 ) 


and is in agreement with reference 1. 

It should be pointed out that the maximum cross-sectional area is in 
general not coincident with, and may be larger than, either of the 
areas A or B (see fig. 1(a)). This result is contrary to two- 
dimensional profiles where the maximum thickness is coincident with the 
largest ordinate through which the profile is required to pass. The 
maximum area of the body of revolution coincides with the point A only 
when 


itA V 1 — C^~ _1/\ / -» r~7 \ 

— = — + COS 1 (-c) (17) 

B c 


The profiles in this case give minimum wave drag for given length, base 
area, and maximum area, and the area distribution is given by 


S(x) = ^ \[l - x 2 + ^ ~~ ~ -loggN + | cos-M-x) (18) 

VI - C 2 


The quantity c in this equation represents the location of the maximum 
diameter and is determined, for any given ratio B/A (A = S max ), from 

equation ( 17 )- A graph of equation ( 17 ) is shown in figure 2 and some 
typical profiles defined by equations ( 17 ) and (18) (which are special 
cases of equation ( 15 )) are presented in figure l(b). 
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The wave drag for the optimum bodies of equation (15) is easily 
found to be 




— ® cos" 1 (-c) 
(l - c 2 ) 2 


Be 

(1 - caH 


+ 


A B 

A - — cos 
n 


-1 


(-C) 


(1 - ca): 


\/l - 4- COS“^(-c) 


B 1 
ix 


+ c \]l - 

(1 - = 2 ) 2 



(19) 


The drag coefficient for optimum bodies of a given length, base 
area, and maximum area is then found to be 




( 20 ) 


where A = S max and c is defined by equation (17). This drag coef- 
ficient is shown in figure 3 plotted against B/A. 

If the optimum profiles which must pass through several specified 
points were sought, it follows from the present section that integral 
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equations like equations ( 11) and (12) must be satisfied in each of the 
specified regions. 

Of some practical interest may be the least-drag body with a given 
length and base area and which must have as a part of its profile a 
cylindrical collar of constant radius. The problem in this case reduces 
to the solution for a tandem biplane configuration having the same upwash 
distribution on the wing surface as in the previous problem. 

Optimum body having given length, volume , and base area .- It follows 
from the previous section that again the optimum bodies having given 
length, volume, and base area will also have a zero slope at the base. 

The calculus of variations then gives the following integral to be made 
stationary: 


h 


-1 -1 


S"(x)S"(!)log e |x - 5| dx d ! + * / S(x) dx 


-1 


( 21 ) 


with the following conditions on the body shape: 


S(-l) = S'(-l) = S » (1) = 0 


S (l) = B 



dx = V 


> 


( 22 ) 


By following a procedure similar to that used in the previous section, 
the following integral equation is deduced: 


S'(x) 
x - I 


dx = K]! 2 


+ Kgl f K 3 


-1 


(-1 < | ^ 1) (23) 
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Again by use of the wing analogy, the solution to equation (23) is found 
in reference 6, and the constants are determined from equation (22). The 
final result is 


S(x) 


8 V - B 
3 * 


(1 



7t V 


B 

- 1 - — cos 

It 


-(-x) 


(24) 


For vanishing base area B this result reduces to the optimum body for 
given length and volume as found in reference 2. 

Typical body profiles as given by equation (24) are shown in fig- 
ure 4. The profile shapes at the nose and base of the body can be seen 
to have essentially the same character as the profiles found in the pre- 
vious section. 

The wave drag for those bodies having a given volume and base area 
is found to be 


_IL_ = l[sv(v _ 2B) +■ 9B 2 ] (25) 

4 ) 


The location of the maximum thickness for the bodies defined by 
equation (24) is given by 



(26) 


Therefore, for a given ratio v/B or a given location of maximum thick- 
ness, the ratio of the base area to the maximum area is readily found. 

A graph of the location of the maximum thickness against Bj S max is 

shown in figure 5- The wave-drag coefficients for the same values of 

B/S are shown in figure 6. 

I max 
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CONCLUDING REMARKS 


By use of an approximate equation for the wave drag of slender 
bodies of revolution in a supersonic flow field, the minim urn -wave -drag 
body shapes have been determined for three cases: (l) the body has given 

length and base area and a contour that passes through a prescribed point 
between the nose and base, (2) the body has given length, base area, and 
maximum area, and (3) the body has given length, base area, and volume. 

The optimum body profiles are shown to have: (l) infinite slope at 
the nose but zero radius of curvature so that the bodies appear to have 
pointed noses, (2) a zero slope at the body base, and ( 3 ) no variation of 
wave drag with Mach number. For those bodies having a specified inter- 
mediate diameter (that is, location and magnitude given) the maximum body 
diameter is found to be, in general, larger than the specified diameter. 
It was also found that, for bodies having a specified maximum diameter, 
the location of the maximum diameter is not arbitrary but is determined 
from the ratio of base diameter to maximum diameter. 

In order to find body profiles that give minimum total drag, con- 
sideration must be given to skin-friction and base-pressure drag. As 
for skin-friction drag, it can be said that the profile shape is of 
secondary importance and that the wetted area of the body and the 
Reynolds number are the primary factors. 

In cases where the base pressure is independent of the body shape, 
for example, where a supersonic jet exits at the base, the optimum pro- 
files sought are the same as those determined herein. If no jet exits 
at the body base, the base pressure is dependent upon the body shape. 

Just what the dependence is is not sufficiently clear at the present 
time. However, the profiles determined may require some slight changes 
and further investigation, in which account is taken for both the wave 
drag and base-pressure drag, would be necessary. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , August 21, 1951 
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(a) Bodies having a given length and base diameter and a contour passing 
through another specifically located diameter. 

Figure 1.- Typical optimum body profiles defined by equation (15). 

Thickness, 20 percent. 
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(b) Bodies having given length, base area, and maximum area (A = Sj^). 

Figure 1.- Concluded. 
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Figure 2.- Location of maximum diameter for optimum bodies of given 
length, base area, and maximum area (A = Sj^), 



Figure 3 .- Wave-drag coefficient for optimum bodies of given length, 
base area, and maximum area (A = S^). 
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B = _1_ 

Smax 16 


e * O.OlU 


C D 

(d /l ) 2 


10.39 




Figure l;.- Typical optimum body profiles which have given length, volume, 
and base diameter defined by equation (21;). Thickness, 20 percent. 
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Figure 5>.- Location of maximum body diameter for a given ratio of base 
area to maximum area for bodies of a given length, volume, and base 
area. 



Figure 6.- Wave-drag coefficient for a given ratio of base area to maximum 
area for bodies of a given length, volume, and base area. 


